Effects of resonant single-particle states on pairing correlations 
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Effects of resonant single-particle (s.p.) states on the pairing correlations are investigated by an 
exact treatment of the pairing Hamiltonian on the Gamow shell model basis. We introduce the 
s.p. states with complex energies into the Richardson equations. The solution shows the property 
that the resonant s.p. states with large widths are less occupied. The importance of many-body 
correlations between bound and resonant prticle pairs is shown. 
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One of current topics is concerned with the microscopic 
structure of nuclei far from the /3-stability line. In such 
nuclei, weak binding of nucleons is expected to cause 
exotic situations where resonant and continuum single- 
particle (s.p.) states contribute considerably to nuclear 
properties even at low energy. The problem has been 
considered since an early time 0] , and recent experimen- 
tal study of drip-line nuclei has stimulated a new inter- 
est in it. The resonant and continuum s.p. states could 
enhance pairing correlations which are most impor- 
tant in nuclei once the Hartree-Fock (HF) or shell-model 
mean field is formed. They could affect binding energies 
of bound states and nuclear resonant states, and could 
play important roles in properties of these low-energy 
states. Recently, theoretical studies 0, |E S 13 
contributions from the resonant and continuum states 
have been reported. Refs. 0, incorporated the ef- 
fect of the resonant continuum or resonant states in the 
HF-I-BCS approximation. Refs. Q studied the effects 
of the resonant and continuum s.p. states on one pair 
state outside a core. 

The purpose of this paper is to investigate the effects 
of the resonant s.p. states on many-body correlations 
in many particle systems (including correlations between 
bound and resonant particle pairs), on the shell-model 
basis. We treat the problem by exactly solving the pair- 
ing Hamiltonian including resonant s.p. states. Namely, 
we extend the Richardson equations ^ so as to include 
s.p. states of complex energies and solve them. The 
Richardson equations have a defect that the pairing in- 
teraction strength cannot be changed for bound and res- 
onant s.p. states. Still, this approach is very useful for 
getting a general perspective of our subject. 

We consider a system of Fermions in which the mean 
field and residual interaction are approximately repre- 
sented by a one-body potential and the pairing force. 
We denote a Fermion by cl^ (cr = ± denote time- 
reversal conjugate states) and a pair operator by Sl — 
c\j^c\_. The eigenstates of our system are classified 
by the number of unpaired Fermions. The states with 
no unpaired Fermion are lower in energy than the oth- 
ers. We pay our attention only to such pair-correlated 
states composed of the basis states • • -jO). In 

this subspace, the pairing Hamiltonian is written as 



H = 2ei,Ni, — GJ2f^i, Sj^S,y, where Si, denotes the s.p. 
energies and G is the pairing force strength. The operator 
Ni, = SlS„ satisfies the relation [N^,Sl] — Sl from the 



Pauli principle cj,g.cj, 







and counts the pair number in 
the level v. Richardson [8| proved that a pair-correlated 
eigenstate is expressed as Ofc^iISi/ '5'J/(2e^ - Zfe)}|0), 
using the roots Zk of the coupled equations 
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Here, the number of the equations ^ is the number of 
pairs N . For an eigenstate, the total energy and occupa- 
tion probability of the s.p. level v are given by 
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Zk, 



k=l 



dG ^ 2ei, - Zk 



(2) 
(3) 



Although Richardson ^] supposed to be real num- 
ber, his proof is not relevant to whether are real or 
complex. We can therefore introduce resonant s.p. states 
with complex energies to our pairing problem by consid- 
ering a Woods-Saxon potential. In compensation for this, 
however, the Hamiltonian stops being Hermitian. Even 
when the s.p. energies Si, are real, some of the roots Zk 
are complex conjugate (the total energy of the system is 
real). If any of are complex in Eqs. (Q), all the roots 
Zk are complex but not complex conjugate. If the res- 
onant s.p. states participate in the pairing correlations, 
therefore, the total energies of the pair-correlated states 
including bound states become complex in general. 

Let us first investigate basic properties of a general 
system which has resonant s.p. states. We consider a 
model that has a dozen s.p. levels with equal spacing 
e, where the lower eight levels are bound states and the 
upper four levels have the same resonant width 7, 



8.5)e for = 1,2,- • -8 , . 

8.5)e- i(7/2)e for = 9, 10, 11, 12. ^ ' 



We use the energy unit e = 1 (also for G and 7) in the 
following pages. 
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When the pair number N is small, the Fermi surface 
lies deeply below the resonant s.p. levels. The resonant 
states can be neglected as long as the pairing force is 
not very strong. This is the case of ordinary treatment 
for the nuclear pairing problem. Calculated results for 
the = 4 system are summarized as follows. The sys- 
tem gets additional energy gain when the resonant lev- 
els are taken into account and the energy gain increases 
more rapidly than that in case of no resonant levels as 
the pairing force becomes stronger. The total energy be- 
comes complex, namely, eigenstates including the ground 
state have finite resonant width when the resonant s.p. 
states contribute to the pairing correlations, as men- 
tioned above. Interestingly, the ground state has pair 
decay width larger than low-energy excited states. 

We are interested in the situation of weak binding of 
nuclei where the Fermi surface lies near the resonant s.p. 
levels. Fig. 1 shows the effects of the resonant s.p. lev- 
els on the = 8 system with 7/2 = 0.5, in which the 
Fermi surface is between the bound and resonant levels. 
If we do not take account of the resonant levels, there is 
no configuration mixing like the closed shell {v^ = 1 for 
J/ < 8) and hence the energy gain due to the pairing cor- 
relations increases monotonously depending on G. Fig. 1 
indicates more rapid increase of the energy gain when the 
resonant s.p. states are taken into account. The imagi- 
nary part of the total energy shown in the lower graph 
of Fig. 1 displays interesting behavior. When G = 0, the 
ground state is the closed-shell state with no decay width 
and the second pair-correlated state is a resonant state 
of the Fermion pair Sl with the half width r/2 = 7 = 1. 
As the pairing correlations become stronger, the "width" 
of the ground state increases and after a critical strength 
of G it exceeds the width of the second state which is 
originated in a resonant state. 

How does the imaginary part of the ground-state en- 



ergy Im(£') depend on the parameter 7? The upper 
graph of Fig. 2 illustrates lm{E) as a function of 7/2, 
in case of iV = 8 and G = 0.3. The value of Im(£^) is 
proportional to 7 in a small 7 region, but reaches the top 
at 7/2 « 2 after becoming a gentle slope around 7/2 ^ 1. 
To see the structure change, we calculated the occupation 
probabilities of respective s.p. states for the ground 
state. The real parts of for the resonant s.p. states 
are illustrated in the lower graph of Fig. 2. This figure 
indicates that the occupation probabilities Re(w^) cross 
the zero point one after another. Especially, Re(w|) of 
the lowest resonant s.p. level v — 9 decreases drasti- 
cally as 7 increases, and becomes negative at7/2R:!l.l5. 
The occupation probability must be positive to keep the 
physical meaning. If one of Re(w^) is negative, the state 
is unphysical. We should, therefore, stop our treatment 
in the unphysical situation for 7/2 ^ 1.15. This model 
calculation suggests that the resonant s.p. states with 
too large widths should be removed from the space of 
the pairing correlations. 
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FIG. 1: Energies of the ground and next pair-correlated states 
of the N = 8 system when 7/2 = 0.5. Variations of real and 
imaginary parts depending on G are illustrated. 



FIG. 2: The imaginary part of the ground-state energy Im(£') 
and the occupation probabilities Re(w^) of the resonant s.p. 
levels as a function of 7/2, in case of A'^ = 8 and G = 0.3. 

It should be noted that Figs. 1 and 2 give information 
about how the quantities vary when the pairing interac- 
tions are weakened for the resonant s.p. states or widths 
7 are narrowed. 

When N — 9, the last pair occupies the resonant s.p. 
levels in our model. This situation is similar to the sys- 
tem of a core plus one pair considered in Refs. 0], 
which discussed important roles of the resonant and con- 
tinuum s.p. states in binding of the last one pair. In 
Fig. 3, we show the total energy of the N = 9 system 
(ground state) as a function of G. The real part Re{E) is 
measured from the base level -64 that is the unperturbed 
core energy E^^]^2(i^ — 8.5). If we regard the N = 8 inner 
part as a core after Refs. Q , the system is reduced 
to one pair moving in the 4 resonant s.p. levels. Its 
energy shown by the dash curve (a) decreases slowly as 
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G increases and becomes negative. The negative energy 
means binding (confinement) of the last pair in the sys- 
tem, according to Ref . . The imaginary part of energy 
Im(£') remains being 1 in the truncated model (a) which 
has only 4 resonant s.p. levels with 7/2 = 0.5. 

If we take the 8 bomid s.p. levels into account, the 
ground-state energy is very much lowered as shown by 
the solid curve (b) in Fig. 3. The energy crosses the 
zero point at much smaller G (~ 0.11) as compared with 
(a). Obviously, the inclusion of the bound s.p. states 
strengthens the correlation energy. This, however, does 
not mean the binding of the last pair, because the 8 pair 
system without the last pair also gets large binding en- 
ergy. The energy difference between the 9 and 8 pair sys- 
tems, which is equal to the pair decay Q-value (we denote 
it by Qpair), is plotted by the dot curve (b') in Fig. 3. 
The curve indicates that the iV = 9 system cannot bind 
the last pair in our model assuming the same potential 
depth for = 8 and = 9. We can say that the correla- 
tions of the last pair in the unbound s.p. states with the 
inner part are very important but whether the last pair is 
bound or not in the system may depend on complicated 
conditions. (It is expected that the 2N -1-1 = 17 system 
is more loosely bound due to the blocking effect than the 
2N = 18 system.) The imaginary part of energy lm{E) 
increases also in the N = 9 system as G increases, in our 
model without the continuum s.p. states. 



TABLE L Occupation probabilities Va of respective spherical 
orbits a in the ground state of ^''O when G = 0.61 MeV. 
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FIG. 3: Ground-state energy of the A'^ = 9 system (7/2 — 
0.5): (a) neglects the bound s.p. states; (b) includes them. 

The previous work Q pointed out important roles 
of the continuum s.p. states in the core plus one pair sys- 
tem. The calculations Q showed that the continuum 
states contribute to the binding of the last pair in the 
system without bound s.p. states like (a). If the contin- 
uum s.p. states are correctly treated, therefore, we can 
expect that the continuum states may contribute to the 
binding of the last pair in the system including bound 
s.p. states like (b) and that the imaginary part of energy 
Im(_E) decreases reasonably. 
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-189.95 


-189.29 


-187.77 


-186.91 



Now let us consider sixteen neutrons of ^''O as an ex- 
ample of realistic system. We neglect the continuum s.p. 
states, because we do not have a reliable treatment of 
them in our method. We generated the single-neutron 
basis for ^^O by a Woods-Saxon potential with the ra- 
dius i?o = 1.2A^/'^ fm, the surface diffuseness a = 0.65 
fm, the potential depth Vb = —58 MeV and the strength 
of the spin-orbit term Uso = 15 MeV 0. The s.p. en- 
ergies obtained are listed in the second column of Table 
I. The energies £a are similar to those obtained with the 
Cwoik code except for the very low energies of the 
2p3/2 and 2pi/2 orbits. The resonant widths are large for 
2p3/2 and 2pi/2- Using the s.p. energies, we solved the 
Richardson equations for the lowest three pair-correlated 
states of ^''O, by varying the pairing force strength G. 
Though we omit figures, calculated result shows the in- 
crease of the pairing correlation energy as G increases, 
which resembles the small G region of Fig. 1. The gen- 
eral features shown by the model Q are confirmed in 

240. 

Table I lists the real parts of the occupation prob- 
abilities and pairing gap A — Tia^a — ^aG{ja + 
l/2)va\/l — v^, for the ground state. Here, we set 
G = 0.61 MeV so that A is nearly the standard value 
12/\/A MeV when 7q = 0. Table I indicates that the 
widths 7a disturb the pairing correlations and pairs hesi- 
tate in jumping to the resonant orbits with large widths. 
The pairing gap A is 2.45 MeV when neglecting 7^, while 
it is 2.09 MeV when we take account of 7^. This feature 
agrees with that of the HF-I-BCS calculation 0|. It is 
important to notice that the resonant orbits with large 
widths participate little in the pairing correlations. This 
property is sharply revealed if the widths of 2^3/2 and 
2pi/2 are larger. When using another parameters for the 
Woods-Saxon potential, we had negative values of Re{vl) 
for 2p3/2 and 2pi/2, which corresponds to the unphysical 
situation pointed out in Fig. 2. 

Kruppa et al. argued that only resonant s.p. states 
with narrow widths should be included, to keep physical 
meaning. If we exclude the orbits 2^3/2 and 2pi/2 fol- 
lowing them, the results become very different as shown 
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in the seventh column of Table I. Excluding 2^3/2 and 
2pi/2 changes the occupation probabihties Re(w^) very 
much. The pairing gap A decreases from 2.09 MeV to 
1.27 MeV. (It is interesting to note that the values of 
Ke{v'^) approach those of the fifth column including 2^3/2 
and 2pi/2, if G is enlarged so as to give the same gap 
A « 2.09 MeV.) The significant reduction means that 
the resonant s.p. states with large widths, if they are in- 
cluded, are little occupied but still contribute to the pair- 
ing correlations. If we neglect the last resonant s.p. orbit 
1/7/2 further, the pairing correlations are much weakened 
as shown in the last column of Table I (A = 0.73 MeV). 
On the ground-state energy tabulated in the bottom line 
of Table I, whether including or excluding the resonant 
s.p. states 2p3/2 and 2pi/2 is more decisive than the effect 
of the resonant widths 7q. Exclusion of 2p3/2 and 2px/2 
reduces the energy gain about 1.5 MeV (from 189.29 to 
187.77), while the resonant widths 7a reduces the en- 
ergy gain 0.66 MeV (from 189.95 to 189.29). The imagi- 
nary part of the ground-state energy lm{E) becomes very 
small -0.09 MeV from -1.12 MeV, if the resonant orbits 



2p3/2 and 2pi/2 with large widths are excluded, which 
complies with the physical demand |^. 

In conclusion, we have investigated the effects of the 
resonant s.p. states on the pairing correlations by an ex- 
act treatment of the pairing Hamiltonian. The present 
study shows that many-body correlations between bound 
and resonant particle pairs are important and that res- 
onant s.p. states with too large widths should not be 
included in the space of the pairing correlations. We 
have not completed the study of the contribution from 
the continuum s.p. states. The interactions related to 
the resonant s.p. states and especially to the continuum 
s.p. states must be appropriately evaluated. Our ap- 
proach using the Richardson equations is not adequate 
to this problem. It is necessary to use different pair- 
ing interactions in the next step. The conclusion that 
resonant s.p. states with very short life time should be 
excluded from the space of the pairing correlations must 
be investigated further in calculations which deal with 
the continuum s.p. states correctly. 
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